Late time evolution of the gravitational wave damping in the early Universe 
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An analytical solution for time evolution of the gravitational wave damping in the early Universe 
due to freely streaming neutrinos is found in the late time regime. The solution is represented by 
a convergent series of spherical Bessel functions of even order and was possible with the help of a 
new compact formula for the convolution of spherical Bessel functions of integer order. 
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I. INTRODUCTION 



Thorough analysis of cosmic microwave background (CMB) radiation provides a unique test of a standard inflation- 
^ ] ary cosmological model [H-Q- While scalar fluctuations of CMB serve as an invaluable source for exploring density 
1—5 i: of matter and radiation and large-scale structure of the universe, f7 H13l . observations of tensor fluctuations of CMB 

' open a window for searching after a signature of gravitational waves |13l - ll6| . 
^f) ,: The CMB observations done by Wilkinson Microwave Anisotropy Probe (WMAP) [l7| generally support theoretical 
predictions based on the standard inflationary cosmological model. The detailed analysis of the experimental data 
' ^ ' , gives more and more accurate values ^18- 20] for the most valuable cosmological parameters such as baryon density, 
qh- total matter density, Hubble constant, and age of the Universe. 

^ ' Independently of WMAP measurements, there is a long quest for a direct observation of cosmological gravitational 
{^JT)' waves Specially designed for this task Laser Interferometer Gravitational Wave Observatory (LIGO) puts a 

I— I ' major effort in this experimental challenge • A direct observation of cosmological gravitational waves would serve 
as a decisive test for validity of the Einstein general theory of relativity in the same way as the Michelson-Morley 
experiment served as a major proof of the Einstein special theory of relativity. 
" As in the experimental case, cosmological tensor fluctuations pose a challenge also from a theoretical side [23| . 
QQ ' Following S. Weinberg we argue that "The particles of both the cold dark matter and baryonic plasma move too 
'[ slowly to contribute any anisotropic inertia. In tensor modes there are no perturbations to densities or streaming 
• velocities, so there are no perturbations to either the cold dark matter or baryonic plasma that need to be followed 
' here." 

. Therefore the only contributions to the anisotropic inertia tensor are due to photons and neutrinos. Further 
CsJ ' simplification comes from the following argument by S. Weinberg (23l - [25j : "The anisotropic inertia tensor is the 
, sum of the contributions from photons and neutrinos, but photons have a short mean free time before the era of 
^ ' recombination, and make only a small contribution to the total energy density afterwards, so their contribution to 
the anisotropic inertia is small. This leaves neutrinos (including antineutrinos) , which have been traveling essentially 
without collisions since the temperature dropped below about 10^" K, and which make up a good fraction of the 
^ ] energy density of the universe until cold dark matter becomes important, at a temperature about 10^ K. The tensor 
- - part of the anisotropic inertia tensor is given by 



n,,{u) = -4p^{u) I K{u - U)K^{U)dU, (1) 

so the gravitational wave equation 







now becomes an integro-differential equation 

j2 



htj + i-^hij - ^htj = IGnCTTij, (2) 



72 • J A ru 

+ 3--rMj - —h,j = -64^Gp,(u) / K{u - U)Kj{U)dUy (3) 
(it d Uit d Jq 
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The impact of neutrino source on gravitational wave damping has been thoroughly considered f23l-l33l| . In this paper 
we report an analytical solution for the damping of gravitational waves in the early Universe due to freely streaming 
neutrinos, eq. (jS]), in the late time regime, u ^ Q ^ 1. The solution is represented by an infinite series of spherical 
Bessel functions of even order. First we shall explain each term in the introduced equations ([T][3]). 



II. NOTATIONS 



We are interested in the time evolution of ft,y (x, that is the tensor perturbation to the metric g^j^^: 

goo = -1, (4) 
9i0 = 0, 

where a{t) is the time-dependent Robertson- Walker scale factor. The kernel K{u) in eq. ([T]) is represented by the 
sum of three spherical Bessel functions j„(u). 

The anisotropic stress tensor 7Tij{u) is obtained from the solution of the Boltzmann equation for freely streaming 
neutrinos [24i [25j which defines the stress-energy tensor for the Einstein field equations. The functions 'p^{u) and 
'Pj{u) give the unperturbed equilibrium neutrino and photon energy density, correspondingly, which define the ratio: 



MO) = = "'''\']'4/3 ^ 0.40523. (6) 

P. + P, 1 + 3(1) (A)^/^ 

The variable u is the product of the wave number k and the conformal time. 
The boundary condition to eq. ([3]) is 

/i^(0)=0, (8) 
and it is assumed that we can parameterize the tensor hij(u) as 

= %(0)x(m). (9) 

Introducing the dimensionless quantity 

y^'-^. (10) 

d^q (l(t — t^q)^ 

where teq is the time of matter-radiation equahty, the general eq. ([3]) can be written as [23 | 

<' - ^) ^ - «^^<'> r - ^)) 

with the boundary conditions 

X(0) = 1, (12) 

dxiy) 



3(1) m 



4/3 



dy 

Eq. (ITT|) can be further simplified by the change of variable. 



= 0. 

v=o 
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into 

Here Q is defined by the ratio of the wave number to its vahie at the time of matter-radiation equality, 

Q^V2^, (15) 

"'eg 

Hit) = ^. 
a{t) 

III. TIME EVOLUTION OF THE GRAVITATIONAL WAVE DAMPING 



In the late time regime, u ^ Q ^ 1, the general eq. simplifies into 



du^ u du Jq ail 

where C = — 24/i,(0)(4(3)^. Below we present the analytical solution for eq. ([T6| together with boundary conditions 

Clearly we need to find a specific function that would "absorb" all derivatives of u and all powers of u on the left 
hand side of eq. (|16p . For the differential operator that appears in the left hand side of eq. p^ . 

f_ 

\ dv? u du 

these conditions can be satisfied with the function 

n(n + 3) ■ ( ^^ , {n-2){n + l) , , . , 
fn[u) = 2n-l + 2n + 3 + ^"+2 WJ- (18) 

Applying the differential operator p7|) to the function ((T8)) we obtain a single spherical Bessel function 

L[fn{u)] - n{n - 2){n + l)[n + 3)(2n + l)j„(u) (19) 

which is exactly what we are looking for. Therefore we can look for the solution of eq. (jl6p in terms of the expansion: 

, , /n(n + 3).. , , . / M (n — 2)(n + l). , , . , ,A 

X{u) = 2]c„ ( l^^_^ bn-2{u)+jn{u)] + 2n + 3 ) ■ (20) 

n— ^ ^ 

The left hand side of eq. (fTG)) transforms into 

oo 

^ n(n - 2)(n + l)(n + 3)(2n + l)c„j„(M). (21) 

The regular at the origin solution for the homogeneous part of eq. (1161) . 

Xo(«) + ^Xo(")+Xo(")=0, (22) 
is the sum of the two spherical Bessel functions, 

Xo{u) io{u) + j2{u). (23) 

The homogenous part Xoi"^) of the general solution x(w) can be already seen as a linear combination of the first two 
terms in the expansion (j20p for n = and n = 2. 

The right hand side of eq. (IT51) is represented by the convolution of the kernel ([S]) with the first derivative of the 
unknown function x(m) which we are looking for in terms of a series (j20p . Clearly one needs a mathematical tool that 
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FIG. 1: Graphical representation for the upper but one triangular structure of the matrix B2k.2i, eq. 
I and k define a position of the matrix element in the ly-plane, while along z-axe we plot its value. 



27l). The matrix indices 



relates a convolution of spherical Bessel functions to a series of those. In Appendix A we prove a useful formula for 
the convolution of spherical Bessel functions, that is not presented in the mathematical literature: 



4(-^)"^ 



2i 



dUj^{u-U)j^{U) 

oc ( 

Y^{2l + l)^ji{u) 



1=0 



^y^"" {l,0,m,0\L,0f 
^ L + L'^-n{l + n) 



1 t=|!-,Tl| 

\ t-n-l>0 



L=l+n 

E 

L=\l-n\ 

7T1 - i - 1 > 



(/,0,n,0|i,0)^ 
L + L"^ — m(l + m) 



where {l,0,m,0\L,0) are the Clebsch-Gordan coefficients. 
The right hand side of eq. ([T6| is 



c r K{u-u)'^du^c f; m B, 



fc,/Cfe. 



L = 



k = 



Here the matrix B^.i is generated by the convolution of the first derivative of the function x(?i), 



n \ 



, , (-2 + n)n(3 + n) , n (-3 - 4n + r?) 



(-3 + 2n)(-l + 2n) 



(-3 + 2n)(3 + 2n) 



. ^ J-1 -n) (2 + 6n + n2) , T-l - n)(-2 + n)(3 + n) 



(-l + 2n)(5 + 2n) 



(3 + 2n)(5 + 2n) 



(24) 



(25) 



(26) 
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2k,2l 



bv 


means 


of eq. 


(p4l) , and for inte 


ger k e [0, 10] and / G [0, 10] is 
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1 
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29 


1 


37 


41 


15 


10 


350 


4900 


3780 


19404 


1981980 


193050 


17381650 


42031990 


1 


1 


13 


17 


1 


5 


29 


1 


37 


41 


3 


2 


70 


980 


756 


19404 


396396 


38610 


3476330 


8406398 


4 


7 


5671 


10523 


733 


1453 


29 


173 


17279 


65887 


15 


55 


17325 


44100 


13860 


640332 


77297220 


5855850 


1095043950 


8322334020 





18 


793 


3013 


13507 


2077 


1247 


353 


1147 


8077 


55 


3465 


8820 


41580 


23716 


220220 


1351350 


56156100 


3627684060 








88 


68323 


6611 


1343 


3652637 


4375681 


5229247 


245303 


225 


209475 


17955 


3234 


30060030 


491891400 


10220410200 


3530687160 











26 


1661 


81703 


4520839 


39881 


67721729 


174227819 


57 


3933 


204930 


8918910 


257400 


5631654600 


231512200920 














12 


74315 


244553 


302161 


62891453 


95878295 


23 


143451 


567567 


504504 


334486152 


6355236888 

















238 
405 


701017 
1142505 


2355893 
5077800 


5773813 
8353800 


7011533 
31744440 




















304 
465 


11539 
16275 


119899 
240975 


1434877 
1831410 























18 
25 


39997 
49725 


1405973 
2645370 


























92 
117 


602003 
669123 



(27) 



/ 

We should notice an unpleasant feature of the matrix of coefficients (P7)) : the first row up to a factor fc = — 5 is 
identical to the second row. This is a direct response to the symmetry of the introduced function /„ («), eq. ^ 
The functions for n — and Ji — 2 are exactly the same up to the factor k = —b, 



f2{u) 
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ijo{u) + j2{u)). 



(28) 



Therefore the rank of the matrix 1^7^ is Rank[i3] = — 1 . This is a real obstacle because it leads to inconsistency 
with the boundary conditions. Indeed, the boundary conditions (|12p are met if we set 



-3C0 



Co = 1. 
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On the other hand, the linear dependence of the matrix is equivalent to 



10 



r^o 



-Co = 0. 



In order to avoid this unpleasant feature we can start summation in the series (I20[) from n 
is equivalent to setting 



Co 



and thus the boundary conditions ()12p are met if we set 



C2 



0, 



3 

To' 



(29) 
(30) 

2 instead of n = which 
(31) 



(32) 



Absence of the j2 (u) term in the left hand side of eq. (|21l) leads to a restriction on the first coefficients in eq. ([^ 



1 

— Co 

3 ^ 



15^^ = 



Finally we get the system of linear equations 

n(n-2)(n + l)(n + 3)(2n- 



C 



n+2 

E 

k = 



(33) 



(34) 
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that returns solution (for even n and k) 

n{n - 2){n + l){n + 3)(2n + l)c„ - C ^" .=o Bk.nCk 

Cn+2 = ■ . (35) 

In the hmit Q ^ 1 and owing to C = — 24/y(0)(4(5)^ we have Q- independent solution: 

k=o Bk^nCk 

Cn+2 = , (36) 



and therefore for n € [0,9] we have 



_ , 3 3 5 35 63 231 429 6435 121551 

- i ' Y^' 8' Ye' 128' 256' 1024' 2048' 32768' 65536 J ^ ' 



which completes our series solution eq. (j20p . 



IV. CONCLUSION 



We have analyzed the problem of gravitational wave damping in the early Universe due to freely streaming neutrinos 
in the late time regime u ^ Q ^ 1. As in the opposite limit u <^ Q [l^, the solution is represented by a convergent 
series of spherical Bessel functions of even order and is independent of the Q-value. Thus we conclude that the 
problem gravitational wave damping in the early Universe due to freely streaming neutrinos is completely solved in 
an analytical way in both early and late time limits. 
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Appendix A 

Here we derive the convolution integral (|24l) of spherical Bessel functions for integer orders n, m. Theoretical efforts 
for such a convolution integral were put forward in (26j but here we report a new compact formula which has a clear 
exchange symmetry {ni — >m) and can be readily applied for a practical calculations. 



Starting with the integral 

Jn,M= I dUj^{u-U)jn{U) (38) 



we prove (1241) . First, we represent a spherical Bessel function as a Fourier transformation of the Legendre polynomial 

jn{u) ^ — 2 — J dsexp{ius)Pn{s). (39) 

Substitution of ^ into leads to 

Jn,mi^)^- — 1 / ds dtexp{iut)Pmis)Pnit) dU exp{iU s - iUt) = (40) 

4 J_i J_i Jq 



I" ds f dtP„(.)P„(i)?"P(^"^)-^"P(^"^) 



4* 

Now we employ the Legendre function of the second kind Qn{z) defined as 



If' , Pjz') 



Performing the integrations over t and s we obtain 



^ ^ ' d5 / dtP^{s)P^{t) / ^ = (42) 



= y diexp(iui) [P„(OQm(i) + /'m(i)Qn(0] • 

Further we reincarnate spherical Bessel functions by decomposing plane waves in terms of Legendre polynomials: 



exp(mt) = ^(2/ + l)^'J^(^i)P^(^), (43) 



(=0 

which leads to 



dtexp{iut) [Pn{t)Qm{t) + P.m{t)Qn(t)] = (44) 



1 



2i 

2i 

1=0 

The angular momentum coupling simplifies the product of Legendre polynomials, 



/ ■\n+rn ^ 

^(2; + l){~l)-'ji{u) / dtPiit) [Pn{t)Qrn{t) + Prn{t)Qn{t)\ • 



L — l-^rn 

Pi{x)P^{x)= J2 {l,0,rn,0\L,0fPL{x), (45) 

L—\l — ni\ 



in terms of the Clebsch-Gordan coefficients 0, m, 0|L, 0). Introducing 



(2/c + l)(m-/c) 



and using the analog of eq. (|45)) , 



Pi(-z)Qm(^)= E (^0,m,0|i,0)'(QL(z) + P^i_i(z))-PK-z)W'™-i(z), (47) 

L = \l-m\ 
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we come to 



dtPiit) [Pn{t)Qrn{t) + Pm{t)Qn{t)] = 



(48) 



dt 



dt 



Pn (t) Pi {t)Q,n (t) + Pi {t)P,n (t) Qn {t) 




Pn{t)\ {l,0.m,0\L,0f{QLit) + WL-iit))-Pi{t)W^^.iit) 

L— |/ — m| 



L—l-\-m 

+ {l,0,m,0\L,0fPL{t)Qn{t) 

L = \l-m\ 

where in the first term we have decomposed the product Pi{t)Qm{t), whereas in the second term we decomposed the 
product Pi{t)Pm{t). Using the parity identity, 



J dtPn{t)QL{t) = - J dtPL{t)Qn{t), 



(49) 



we obtain 



dtPiit) [Pn{t)Qra{t) + Pra{t)QrM 



(50) 



dt 



L—l-\-m 



Pn{t)\ Y {l,0,^n,0\L,OfWL-l{t)~Pl{t)W^^i{t) 
L=\l-m\ 



Therefore 

L—l+m 

dt Y {l,0,rn,0\L,OfPr,{t)WL-i{t) 



(51) 



L^ll-ml 
^1 L— /+m 



[(L-l)/2] 



dt Y (^0,m,0|i,0)2 Y 

"1 L=\l-m\ k=0 



L—\l — 771 

On the other hand, 



2i - 4fc - 1 
(2fc+ 1)(L- fc) 



L—l-\-Tn 



PL-2k-l{t)Pnit) ^ 4: Y 



i-Tl-l>0 



(;,0,m,0|L,0)^ 
L + L'^ - n{l + n) ' 



dtPi[t)W„,-l{t)Pn{t) 



(52) 



[(m-l)/2] 



//' s 

•'^^ k=0 
1 [(™-l)/2] 



2m - 4fc - 1 
(2/c + l)(m- fc) 

2m - 4fc - 1 



Prn-2k~lit)Plit)Pnit) 
L=l+n 



fc=0 



L=\l-n\ 



m-L-l>0 



L + — m(l + m) 



Thus, the identity ([M)) is demonstrated 

^1 L — l+m 



/I L — l+m „1 
rff E (Z,0,m,0|L,0)'P„(t)M^L-i(t)- / dtP„(t)P;(i)W^m-iW = 
-1 r_l; ™l "'-I 



(53) 



^ (/,0,m,0|L,0)' ^ ^ (/,0,n,0|L,0)' 



L=|i-m| 



7Tl-i-l>0 
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